Dynamical Response of Single Bi-layer Spin Model : A Theoretical Analysis 
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The spin dynamics of a single bi-layer ferromagnetic model, as proposed for a manganite system 
La\.<2 1 Sr\&Mn<2 1 0*i as an approximate minimal model, is studied both below and above the bulk 
magnetic transition temperature using the semi-classical Monte Carlo-molecular dynamics tech- 
nique. The quantities studied are:- (i) the static spin configurations, (ii) the spin auto-correlation 
function (SAF) C(t) and (hi) the dynamical structure function (DSF) S^q, uj). The major aim is to 
probe the nature of collective modes and in particular to ascertain whether the propagating modes 
can exist above the ordering temperature. We find that the typical spin configurations contain a 
high degree of mis-aligned spins, particularly in the high temperature phase. Nevertheless, in the 
low temperature phase a long range ferromagnetic ordering can be seen. The shape of the curves 
corresponding to S^q, uj) vs. uj in the constant q-scan are found in general to be quite complex viz. 
containing peaks, cusps and broad plateau. Moreover, in a large regime of (q, uj) space, S^q, uj) 
is found to be negative, signalling a total breakdown of the semi-classical approximation in the 
presence of the enormous thermodynamic and quantum fluctuations. Considering the physically 
allowed regime only, we can theoretically extract an effective uj vs. q dispersion curve for the col- 
lective excitations, which exhibits a slope opposite to that expected for a full three-dimensional 
system and disagrees with the results of inelastic neutron scattering experiments in the nature of 
the dispersion curve along the < 100 > direction in the ordered phase. The theoretical results also 
exhibit the existence of propagating modes even in the high temperature (disordered) phase from 
our calculations based on a minimal single bi-layer ferromagnetic model. 

PACS numbers: 61.05.fm 65.40.De 



I. INTRODUCTION 

The problem of understanding the spin dynamics of 
the spin models, both classical and quantum, at temper- 
atures above the ordering temperature is a long stand- 
ing challenging problem in condensed matter physics [1- 
4]. There are key questions regarding the nature of the 
collective modes viz. whether they are propagating or 
damped. This problem has assumed greater significance 
and importance in the context of the recent experimental 
studies on the various types of low-dimensional and lay- 
ered magnetic systems like Manganites, Cuprates, Bro- 
mides etc [5-7]. All of these systems are modelled by 
either ferromagnetic or antiferromagnetic type of Heisen- 
berg spin Hamiltonians with the presence of either Ising 
or XY-like anisotropy in a few cases. Furthermore, most 
of these systems have both intra-layer and comparatively 
weaker inter-layer exchange interactions and this makes 
the study of spin dynamics even more challenging. 

In this short communication, we present our calcu- 
lational results for the spin dynamics with parameters 
appropriate to those of Lai^Sri^M^Oj , both above 
and below the bulk ferromagnetic ordering temperature. 
Detailed input from crystallographic and other measure- 
ments leads to the conclusion that in this system the 
intra-bilayer exchange interaction is much much larger 



than the inter-bilayer interaction. This motivates us 
to effectively model the magnetic properties of the sys- 
tem by a single bi-layer spin Hamiltonian. Furthermore, 
the intra-bilayer interaction has two components viz. the 
intra-layer interaction represented by J\ and the inter- 
layer one represented by J2, where the magnitude of J\ 
is larger than that of J2 and both are ferromagnetic in 
nature. We carry out detailed Monte Carlo-molecular 
dynamics (MCMD) studies on the above model by con- 
sidering a lattice of finite size with dimensions 50 x 50 x 2. 
This size then represents a single bi-layer system with size 
of each layer being 50 x 50. The calculations were per- 
formed in the low temperature phase (16K) as well as 
in the high temperature phase (200K). The experimen- 
tal value for the Curie temperature (T c ) of this system is 
128K [3,4]. Various static and dynamic properties were 
studied. The important amongst them are :- (i) the static 
real space spin configurations, (ii) the dynamic structure 
function and (iii) the spin auto-correlation function. We 
try to present a detailed microscopic analysis of these re- 
sults and compare them with the available experimental 
results. The experimental results available so far, are all 
based on the inelastic neutron scattering measurements 
carried out on the colossal magnetoresistive bilayer man- 
ganite Lai^Sri&M^Oj . These have direct bearing on 
our theoretical results for the dynamic structure func- 
tion S^q, uj). It would however be equally important 
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and interesting to compare our results for the spin auto- 
correlation function with the corresponding results from 
Perturbed Angular Correlation (PAC) or ESR or Muon 
Spin Resonance experiments for the above system, when 
available. 

The theoretical study of static spin configurations is 
also extremely challenging and gives us a lot of useful in- 
formation regarding the intrinsic fluctuations originating 
from the spin model. These again in general, can have 
contributions from the quantum effects, the thermody- 
namic effects and also spatial dimensionality as well as 
the finite size effects. It is however very diffficult to ex- 
perimentally determine the static spin configurations in 
real space directly. Nevertheless, the real space spin con- 
figurations do play a very crucial role in determining the 
nature of both the static and dynamic spin correlation 
functions. 



II. MODEL AND BRIEF CALCULATIONAL 
PROCEDURE 

The effective spin model we employ is essen- 
tially a nearest and next nearest neighbour quasi- 
two-dimensional (describing a single bi-layer) quantum 
Heisenberg spin Hamiltonian, given by [3] 

T-L = — [Yl < ij > J\Si ■ Sj + T l<< ij >> J2Si - Sj] (1) 

where J\ and J2 are the nearest and the next nearest 
neighbour ferromagnetic exchange interactions respec- 
tively. Ji operates within a layer whereas J2 operates 
between the two layers of the bilayer. The above model 
corresponds to a single bi-layer, as explained earlier. It 
may further be mentioned that the magnetism arises from 
the electrons of Mn atoms possessing an effective spin 
S = 1.8 [3]. The values of J\ and J2 are 4.8 mev and 1.7 
mev respectively [3,4]. 

We carry out a semi-classical MCMD calculations on 
the above model. The details of this technique can be 
seen in the earlier works [1]. The semiclassical treatment 
is expected to hold for the present model, as the effective 
spin value per site is close to 2. 

We implemented Metropolis algorithm for cooling, 
starting from infinite temperature. Corresponding to 
T = 16if , we stored spin configurations with MC ages 
between 5000 and 15,000 MCS (MC steps per spin). For 
200K however, we stored spin configurations with MC 
ages upto 30,000 MCS. 

We carried out the tests for stabilization of various 
thermodynamic quantities to ensure the attainment of 
thermal equilibrium. To our surprise however, we dis- 
covered that in our present model system even when the 
above stabilization or convergence sets in, the fluctua- 
tions in the static properties are unexpectedly large. We 
address this issue in the section on Results. 



III. RESULTS 

In this section, we present our calculational results for 
both static and the dynamic properties. Moreover we 
highlight the salient features of our results. 

Let us first discuss about the static properties. We 
analysed the real space spin configurations and calcu- 
lated the average magnetization, average energy and the 
specific heat from these configurations. This was done 
for both the temperatures, T = 16K and 200K. 

We observed that in the low temperature phase (16K) 
the most probable spin configurations generated by MC, 
have a moderate degree of misalignment of spins, viz. 17 
percent. The average magnetization per site is however 
almost constant. The ratio of the square root of the 
specific heat to the absolute value of the average energy 
is close to 0.4, which is much larger than the expected 
magnitude for our system in the thermal equilibrium. 
This brings into focus the substantial and non-trivial 
contribution of the fluctuations to the static properties 
even in the low temperature phase. 

In the high temperature regime viz. at 200K, our MC 
data show that the misalignment of spins in a typical 
spin configuration is around 80 to 90 percent. Moreover 
the average magnetization at a site also exhibits a mis- 
alignment (disorder) of degree 40-60 percent. The more 
astonishing result is found in the heat capacity calcula- 
tions. The ratio of the square root of specific heat and the 
absolute value of the average energy attains a magnitude 
between 2.0 and 3.0 whereas the expected magnitude is 
around 0.02 in thermal equilibrium at any temperature 
for our system size. 

The above analysis involving the static properties from 
MC calculations brings out the presence of enormous fluc- 
tuations even in the thermal equilibrium. We would now 
like to present and analyse the results for spin dynamics, 
as obtained from our MCMD calculations and examine 
the consequences of fluctuations further. 

The results for the dynamical structure function and 
the spin auto-correlation function show quite different 
type of behaviour with regard to the effect of fluctuations, 
both in the low temperature and in the high temperature 
regime. In particular, the results for 5(q, uj) exhibit very 
peculiar and dramatic nature. In contrast, the results 
for C(t) are more well behaved and do not display any 
surprising feature. 

The most interesting and intriguing aspect of the re- 
sults for spin dynamics is the appearance of the negative 
values for 5(q, uj) in a very large region of the u;-space. 
The details of these striking features are given below. 

At the temperature corresponding to 16K, in the 
5(q, uj) vs. uj plot for q along the < 100 > direction the 
positive regime dominates over the negative one, partic- 
ularly in the low q and low uj regime (see fig. [I] (top 
panel)). In the physically admissible regime viz. the 
regime with positive values of cj, there are signatures of 
propagating modes (see fig. [TJ middle panel). The prop- 
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FIG. 1: (Top and middle panel) 5(q, uj) vs. uj for q = 0.1 rlu 
and q = 0.5 rlu along < 100 > direction at T = 16K; (bottom 
panel) C(t) vs. t at T = 16K. 



FIG. 2: (Top and middle panel) #(q, uj) vs. uj for q = 0.1 
rlu and q = 0.5 rlu along < 100 > direction at T — 200K; 
(bottom panel) C(t) vs. t at T = 200K. 



agating modes are characterized by the occcurrences of 
peaks at finite uj in the above plots. The quantity C(t) 
shows a monotonic decay with time (see fig. [TJbottom 
panel) ) . 

In the high temperature regime i.e. at 200K, we ob- 
serve a complete reversed trend in 5(q, uj). The re- 
gion with negative values for DSF dominates over that 
with positive values for almost all the g-values along the 
< 100 > direction ! Therefore it is conceptually dif- 
ficult to determine the nature of the collective modes 
whether they are propagating or diffusive. Furthermore, 
the shapes of the curves representing 5(q, uj) vs. uj in the 
constant g-scans display a lot of variety with the occur- 
rences of cusps, plateau and shoulders besides peaks (see 
fig. [2ftop panel) and fig. ^ middle panel)). The SAF 
shows a monotonic but very sharp decay with time (see 



fig. [2|bottom panel)). The high temperature DSF vs. uj 
plot in the entire q space along the < 100 > direction, 
can also be seen in the three-dimensional color graph (see 

The above two temperature regimes may be character- 
ized as "moderate fluctuation driven" (for the low tem- 
perature phase) and "strong fluctuation driven" (for the 
high temperature phase), on the basis of the analysis of 
our static results. Our spin dynamics results reaffirm 
and manifest these effects more clearly. The "unphysical 
region" observed in the calculations for DSF originates, 
we believe, from the breakdown of the semi-classical ap- 
proximation (used in our MCMD) in the presence of the 
unusually high degree of fluctuations, as we discuss be- 
low. 

It is well known in literature [8] that 5qm(q 5 ^) (DSF 
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as obtained quantum mechanically) is real and positive 
whereas Scl(^ w) (DSF as obtained classically) can 
in general be a complex quantity. Furthermore, the 
dynamical spin fluctuations and the imaginary part of 
the dynamical spin susceptibility, as calculated quantum 
mechanically, can be directly related by the "fluctuation- 
dissipation" theorem. The dynamical spin fluctuation 
reduces to the dynamical structure function (DSF) in 
the paramagnetic phase. Besides, DSF for positive lj 
and negative lj are related by the "detailed balance" 
condition in thermal equilibrium. These features and 
properties prompted Windsor [1,9] to propose an ansatz 
by which one can extract the true quantum mechanical 
DSF from the one which had been obtained using the 
nearly classical methodology of MCMD. It should be 
kept in mind that this ansatz is ideally suited for the 
paramagnetic phase ; however it may be used even in 
the ordered phase provided the magnitude of the order 
parameter is not too large. 

Mathematically, the ansatz of Windsor is expressed as 



^QM(q,w) 



2RealScL(q., 



(2) 



The above prescription takes into account all the above 
mentioned properties and is essentially based on an 
ansatz that the classical spin-spin correlation function 
is equal to the real part of the corresponding quantum 
function in the time domain. This assumption is mean- 
ingful and accurate in the semi-classical limit i.e. gen- 
uinely valid in the high S and in the rather low energy 
(low lj) limit. In fact, it can be shown that as long as the 
energy under consideration is less than the characteristic 
energy (t^lj c ) corresponding to the "natural time unit" 
(to) occurring in our formalism, Windsor's ansatz works 
well. It may be recalled that (t^ 1 ) is proportional to 
J[S(S+ 1)] 2 and this has dimension of energy. Thus it is 
obvious that this parameter lj c increases with S, implying 
a very wide regime for the validity of the semi-classical 
ansatz in the large spin limit. The microscopic origin of 
the existence of this characteristic/cutoff energy scale is 
the genuine semiclassical-like behaviour of the spin dy- 
namics, prevailing in the time scale above the quantum 
fluctuation (spin flip) time to. 

The manifestation of the above property can be seen 
in our MCMD results on spin models with parameters 
corresponding to the Rare Earth Chalcogenides like EuO 
and EuS. The cutoff energies turn out to be 4.5 meV 
and 0.84 meV for EuO and EuS respectively. Above 
these energies the respective DSF's calculated by our 
semiclassical MCMD for the above two systems, when 
plotted in the constant g-scans as a function of lj, indeed 
become negative as can be seen from the corresponding 
figures [1,2]. Thus we may say that these cutoff energies 
are highly physical and truly represent the boundary 
between the semiclassical and the quantum behaviour 
of the spin dynamics. Furthermore, our theoretical 
estimates of these energies based on the arguments 



x10 D 




20 



0.1 40 



Energy transfer (meV) 



FIG. 3: Three-dimensional colored plot of high temperature 
5(q, uj) vs. uj for various values of q along < 100 > direction. 



presented in the earlier paragraph, are also fairly close 
to the observed ones from the plots [1,2]. 

Now coming to the case of our present system viz. 
Lai^Sri.sM^Or, we extract the values of the parame- 
ter lo c from the S(q, lj) vs. lj plots obtained from MCMD 
and also estimate its magnitude from the value of to. 
When we compare the two, we see a very large discrep- 
ancy! Keeping in mind that to = r, we get 

J F & 2ttJi [5(5+1)] i ' 

the magnitude of to as 6 x 10 -14 seconds. This leads to 
the magnitude of the critical energy as high as 59 meV 
approximately ! Our MCMD results however indicate a 
negative regime for DSF even below 5 meV at both 16K 
and at 200K. 



In order to understand the above peculiarity, we will 
have to also include the possible effects arising from the 
thermal fluctuations, seen earlier in the studies of the 
static properties from MC calculations. The giant fluc- 
tuations observed in the thermodynamic properties imply 
that a wide variance is also expected in the contributions 
to the various static and dynamic response functions from 
the different ensemble configurations. The quasi-two di- 
mensionality of the system is the root cause of these gi- 
ant thermodynamic fluctuations. Thus the total fluctua- 
tions in the system far exceed the amount expected from 
the pure quantum effects. Furthermore, the existence of 
these unusually large thermodynamic fluctuations even 
away from the transition temperature may also cause a 
departure from the standard fluctuation-dissipation re- 
lation, leading to the possible modification of Windsor's 
Ansatz. As a result, the extracted Sqm(*1,v) from the 
semiclassical MCMD results (by using the above ansatz), 
can exhibit a lot of unphysical behaviour and may even 
become negative ! 
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FIG. 4: (Blue circles) Dispersion cj q vs. q along < 100 > 
direction for the low temperature phase viz. at T = 16K. 
(Red circles) Dispersion cj q vs. q along < 100 > direction for 
the high temperature phase viz. at T = 200K. 



IV. CONCLUSIONS 

In this short communication we have tried to present 
an analysis and a possible interpretation of our cal- 
culational results for 5(q, uS) and C(0,£). In particu- 
lar, we have put forward an explanation for the occur- 
rences of negative values for the extracted Sqm{^ 
in a large regime of (q, uj) space. The calculated spin 
auto-correlation function from our semiclassical MCMD 
on the contrary, behaves in an expected way — exhibiting 
a faster decay with time at higher temperature. This con- 
trasting behaviour between DSF and SAF arises because 
the former is a "global property" whereas the latter is a 
"local property" and that the presence of a high degree 
of misalignment of spins (seen in the giant fluctuations) 
affect the former much more than the latter. 

Furthermore, our quantum mechanical DSF is calcu- 
lated by using an ansatz which may not work very well 



when the thermodynamics is dominated by the very large 
fluctuations, as is the case here. On the other hand, our 
SAF is computed directly from the MCMD results and 
is of semiclassical nature. Therefore, our calculated and 
extracted estimates of the above two quantities do not 
obey the theoretically expected consistency relation be- 
tween them. 

In the physically admissible regime, the 5(q, uo) vs. uo 
plots can be used for studying the nature of excitations 
and understanding the genuine spin dynamics. The dis- 
persion curves cj q vs. q corresponding to the "propagat- 
ing regime" from our calculations in both the low tem- 
perature phase and in the high temperature phase have 
been presented in the figure [4] These may also be com- 
pared with the corresponding experimental results, for 
the judgement on the rationale of our modelling itself. It 
is quite striking that the slope of the calculated disper- 
sion curve is negative, i.e. opposite to that expected for 
a full 3D system [3]. In order to understand this, it is 
to be noted that the long range order is quite "fragile" 
in this case and hence the collective propagating modes 
can show quite different behaviour in its dispersion from 
that in a pure 3D system. In particular, the propagating 
modes become more ill-defined with increase in the mag- 
nitude of "q" because of "disorder" present in the short 
length scales within the spin configurations in the or- 
dered phase, causing a softening of the modes as well. In 
addition, the inadequacy of the semi-classical treatment 
can cause more deviation from the expected behaviour 
in S(q,uj) results, as obtained by MCMD calculations. 
These lead to the drastically different slope of the dis- 
persion curve of the propagating modes in the ordered 
phase in our MCMD calculations, from that for a full 
3D system. In the paramagnetic phase, the behaviour is 
even more complicated as there is a competition between 
the "short range ordering" expected for a full 3D system 
and "fragility" inherent within the spin configurations in 
a quasi-2D system. This causes the dispersion curve to 
first have a rise and then a fall. 

Before concluding, it should be emphasized that 
the presence of this high degree of thermodynamic 
fluctuations in the static properties seen from the 
MC calculations after the convergence, is not at all 
inconsistent with the attainment of "steady state" for 
the spin system. This state however does not always 
ensure complete "thermal equilibrium" for the system, 
as is understood in the conventional way. Nevertheless, 
even in the complete thermal equilibrium condition 
as is well known, the magnitude of fluctuations can 
very well exceed the usual expected estimate viz. -^tf 
under certain conditions, like in the vicinity of phase 
transition for example. In our case however, the low 
dimensionality viz. the quasi- two dimensional nature 
of the system (lattice) appears to be the most likely 
origin of these enormous amount of thermodynamic 
fluctuations occurring both in the ordered as well as in 
the disordered phase. 
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